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ABSTRACT 

We study implications for the apparent alignment of the spin axes, proper-motion directions, and 
polarization vectors of the Crab and Vela pulsars. The spin axes are deduced from recent Chandra 
\ X-ray Observatory images that reveal jets and nebular structure having definite symmetry axes. The 

alignments indicate these pulsars were born either in isolation or with negligible velocity contributions 
from binary motions. We examine the effects of rotation and the conditions under which spin-kick 
alignment is produced for theoretical models of neutron star kicks. If the kick is generated promptly 
O I during the formation of the neutron star by asymmetric mass ejection and/or neutrino emission, then 

the ahgnment requires that the protoneutron star possess, by virtue of the precollapse stellar core's spin, 
an original spin with period Ps much less than the kick timescale Tkick, thus spin- averaging the kick 
forces on the star. The kick timescale ranges from 100 ms to 10 seconds depending on whether the 
kick is hydrodynamically driven or neutrino-magnetic field driven. For hydrodynamical models, spin- 
■ kick alignment further requires the rotation period of an asymmetry pattern at the radius near shock 

' breakout 100 km) to be much less than Tkick ^ 100 ms; this is difficult to satisfy unless rotation plays 

O^l , a dynamically important role in the core collapse and explosion (corresponding to Pg 1 ms). Aligned 

kick and spin vectors are inherent to the slow process of asymmetric electromagnetic radiation from an 
\ off-centered magnetic dipole. We reassess the viability of this electromagnetic rocket effect, correcting a 

factor of 4 error in Harrison and Tademaru's calculation that increases the size of the effect. To produce 
a kick velocity of order a few hundred km s~^ requires that the neutron star be born with Pg ~ 1 ms 
\ and that spindown due to r-mode driven gravitational radiation be inefficient compared to standard 

magnetic braking. The electromagnetic rocket operates on a timescale of order 0.3 (P/10^'^ G)~^ year. 
The apparent spin-kick alignment in the Crab and Vela pulsars places important new constraints on 
'p^, each of the mechanisms of neutron star kicks we consider. 

O ' Subject headings: stars: neutron - supernovae: general ~ pulsars: individual (PSR B0531+21, PSR 

|H \ B0833-45) - stars: rotation 

1. INTRODUCTION lutionary studies of neutron star binary population also 
: It has long been recognized that neutron stars have existence of pulsar kicks (e.g., Dewey & Cordes 

X ; space velocities much greater than those of their pro- 7; Fryer & Kalogera 1997; Fryer et al. 1998). Fmally, 

, n i i J- r 1 i- • there are many direct observations oi nearby supernovae 

?H genitors. Recent studies ot pulsar proper motion give , / , ,„„„s , , i • i 

rrt ' onn cnn i -1 j-u oT-i i -j. f / (e.g., Wang et al. 1999) and supernova remnants which 

. >T 1 200 — 500 km s as the mean 3D velocity ot neutron , , , • ; i • n 

, J. u • J.1 c T IT- 1 nr> J T • show that supernova explosions arc not spherically sym- 

stars at birth (e.g., Lyne and Lorimer 1994; Lorimer et 'c 

et al. 1997; Hansen & Phinney 1997; Cordes & Cher- . , „ i , • , • • i ;i 

cc inno\ J ■ -c i- i -i-- u • i • While the evidence lor such kicks is unequivocal, the 

not: 1998), and a signihcant population having veloci- ..... . , , , , , . , , , , 

,. , 1 nnn 1 -1 t-v- i -J r 1 physical Origin remains unclear. A natal kick could be 

ties greater than 1000 km s ^. Direct evidence tor pul- ,,i ,• 

1 ■,. innn i -1 r u j.- generated by an asymmetric cxplosion duc to global hydro- 

sar velocities ^ 1000 km s comes from observations ° -i ,, ;- •,, /^i, -, 

r ,1 1 11 J J 1, ii /-I -J. AT u 1 1 dynamical perturbations m the supernova core CrOldreich, 

ot the bow shock produced by the Guitar Nebula pul- t • o n i ^■ ^r^r^n n o tt i^,^^ t • t • 

sar fB2224+65) in the interstellar medium fCordes Ro- ^ Sahrlmg 1996; Burrows & Hayes 1996; Lai 1999; Lai 

DCLl I J_J jLi ZjjLi^ I \JiJ I 111 LllO lllLi'^l O LVjlldrl lllL^Lll Lllll I V_.'Wl LlOO • 1. VjU r. ^ i i • ^ c\r\r\r\ 1 \ • • ill iir i' 

• n T 1 inno\ A 4- 1 1 1- r i. & Coldreich 2000a, b). Or it could be a result 01 asymmetric 
mam & Lundgren 1993). A natural explanation tor such , . . . .' v, „ , 

1 . 1 1 ., . . ,1 , , . , neutrino emission m the presence ot superstrong magnetic 

high velocities IS that supernova explosions are asymmet- nll/n^1^l^^^•.l ; . . / x-o 

• J -J 1-14. 4. 4. 4. T^i 1 faelds (B ^ 10 Cr) m the proto-neutron star (e.g., Lai & 
ric, and provide kicks to nascent neutron stars. Ihegeode- „. , » c t • i i <• ,i ■ \ a 
,. ■ • ii, u- 1 nco inio,ic?i^r- Qian 1998; Arras & Lai 1999a, b and references therein). A 
tic precession m the binary pulsar PSR 1913-1-16 (Weis- ^ , . ^ ^ ■ i • , i- i i tt • i i 
1 T3 • p rp 1 1 non J 4-11 nnn t/ post-natal kick mechamsm studied by Harrison and lade- 
berg, Romani & laylor 1989; Cordes et al. 1990; Kramer ^ i- , • i , i- 
inno TIT 4- 1 r,nnn\ 4-T, u-4- i i • • maru (1975) relies on asymmetric electromagnetic radia- 
1998; Wex et al. 2000), the orbital plane precession m ,. ' ^ , i i- i • -n , ,• ^ 
PSR J0045-7319/B-star binary and its fast orbital decay ^^^^^^^ off-centered dipole m a rapidly rotating pul- 
/TA • 4. 1 innc T • 4- 1 innc r • innG\ j 4-1, i,- i, sar. All these possibilities have intrinsic uncertainties (see 
(Kaspi et al. 1996; Lai et al. 1995; Lai 1996), and the high „„ i n , i i % t^ ; r .i ,i , •, • i-rr. , 

4- • v c TD /V ■ ^ J Tj 1 s3 and q4 below). Part of the reason that it is difficult to 

eccentricities ot Be/X-ray binaries (see van den Heuvel % r , . , . . , i • i i i r 

c r, 1 aa-tn 11 i ii 4-- r i determine which mechanism is at work is the lack ot any 

& van Paradiis 1997) all support the notion ot neutron . , , ,. i • i , i .; i ;i ,. 

, 1 . 1 rpv. ii J 4. 4. J u- u 4- • IV simple relationship between velocity and other properties 

star kicks. Ihe recently detected high systemic velocity ^ ^ , , , t^ i i •, 

/ 1 -i\ r 4.1 V u- r^- ■ V 1 /T u 01 neutron stars (but see q2). boi example, despite some 

(~ 430 km s ^) ot the X-ray binary Circmus X-1 (John- , , . ,,.,.11,1 . ■ 1 r 

. . 1 innn\ iu j jt, 4- 4- early claims, statistically there is no evidence tor a cor- 

ston et al. 1999) can only be produced by a neutron star 1 1 ; 1 •, 1 ;• , r i- 

kick of at least ~ 500 km s'^ (Tauris et al. 1999). Evo- ""^l^*^"^ b^*^^^^ ^^1°^^*^ ^'^^ magnetic moment for radio 
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pulsars (Lorimcr ct al. 1997; Cordcs & Chcrnoff 1998), 
or between velocity and rotation (Dcshpandc, Raniachan- 
dran & Radhakrishnan 1999). Unfortunately, given the 
large systematic uncertainties, these statistical results, by 
themselves, can not be used to support or rule out partic- 
ular kick mechanisms (see §2.3). 

Recent observations of the Vela pulsar and the surround- 
ing compact X-ray nebula with the Chandra X-ray Obser- 
vatory reveal a two sided asymmetric jet at a position angle 
coinciding with the position angle of the pulsar's proper 
motion (Pavlov et al. 2000). The symmetric morphology 
of the nebula with respect to the jet direction strongly sug- 
gests that the jet is along the pulsar's spin axis. In §2.1 
we examine the polarization angle of Vela's radio emission 
which corroborates this interpretation. Evidence for spin- 
velocity alignment also exists for the Crab pulsar (§2.2). 
We argue in §2.3 that current observations of other pulsars 
neither support nor rule out any spin-kick correlation. 

The alignment between the projected spin axis and 
proper motion (for Vela and Crab) immediately demon- 
strates that binary break-up along with symmetric super- 
nova explosions (see Iben & Tutukov 1996) is not likely 
to be responsible for the observed proper motion. Binary 
disruption would yield orthogonal spin and velocity vec- 
tors if pre-explosion binaries have aligned spin and orbital 
angular momenta (e.g., by virtue of tidal coupling) and if 
the neutron star's spin is in the same direction as the pro- 
genitor's. Projection onto the plane of the sky could fortu- 
itously yield apparently parallel spin and velocity vectors, 
of course, but this is a low probability occurrence. The 
observed spin- velocity alignment in Vela and Crab pulsars 
then indicates that the pulsar was born in isolation, or the 
progenitor binary was so wide that the orbital speed was 
negligible compared to the kick. In §3 and §4 below, we ig- 
nore the binary possibility, and examine how rotation may 
affect kick and under what conditions spin-kick alignment 
is expected for different kick mechanisms. 

2. OBSERVATIONAL CONSTRAINTS ON THE GEOMETRY 
OF PULSAR SPIN AND VELOCITY 

Here we briefly discuss observational results related to 
spin - proper motion correlation in pulsars. 

2.1. Vela Pulsar 

The proper motion of the Vela pulsar (/i^, ^5) = (— 48± 
2,35 ± 1) mas yr"-'^ (Bailes et al. 1990) is at position 
angle V/* = —53.9° ± 1.4° and implies a 2D velocity 
~ 70 - 141 km s'^ for a distance D = 0.25 - 0.5 kpc. 
The proximity of the pulsar to the Earth implies that any 
correction for differential galactic rotation is small. The 
X-ray jet and bow-shock axes of the pulsar nebula lie in 
the same direction as the proper motion to within 5 de- 
grees (Pavlov et al. 2000). If the jet originates from the 
pulsar magnetosphere, as seems likely, it is most natural 
to associate the jet axis with the pulsar spin axis. 

Another constraint on the geometry comes from radio 
polarization. The Vela pulsar formed the basis for the ro- 
tating vector model (Radhakrishnan & Cooke 1969), which 
relates the polarization position angle to the rotating dipo- 
lar magnetic field of the neutron star. At the pulse cen- 
troid, the spin axis, magnetic moment and line of sight 
are coplanar. Because the radiation is highly beamed, the 
polarization angle at the pulse center is determined by the 



orientations of the stellar magnetic field and of the pro- 
jected spin axis. Phenomenologically, radio emission often 
conforms to this picture, but the polarization angle can be 
in one of two orientations separated by 90 degrees. Among 
various pulsars, the two modes vary stochastically from 
pulse to pulse as well as systematically across pulse phase. 
In some objects, the emission seems to be dominated by 
one mode in all pulses and at all pulse phases. 

Using data from Deshpande et al. (1999), which have 
been corrected for Faraday rotation in the interstellar 
medium and in the ionosphere, the pulse centroid po- 
larization angle is V'poi = 35° ± 10°. The difference an- 
gle between the proper motion and polarization angle is 
Alp = ippoi — ipfi = 89° ± 11°, i.e., the polarization vector 
is perpendicular to the proper motion axis. 

We conclude that the space velocity and X-ray jet axes 
are parallel in projection on the plane of the sky. Moreover, 
the radio polarization angle is 90 degrees with respect to 
this axis, implying that the polarization mode dominant 
in the Vela pulsar is one where the electric field in radio 
emission is orthogonal to the magnetospheric field. 

The a priori chance of alignment (or anti-alignment) of 
three direction vectors (corresponding to proper motion, 
perpendicular polarization, and X-ray jet axes) is small. 
We assume each is uniformly distributed on a sphere and 
use measurement uncertainties of 1.4, 10, and 5 degrees, 
respectively to infer a probability of 0.6%. The corre- 
sponding probability increases to 1.2% if alignment of ci- 
ther the parallel or the perpendicular polarizations are 
considered. 

2.2. Crab Pulsar 

Caravco & Mignani (1999) report a new HST-derived 
proper motion for the Crab pulsar, {iJ,a,Hs) = (— 17 ± 
3, 7 ± 3) mas yr~^ , or /z = 18 ± 3 mas yr~^ with a position 
angle V'/^ = 292 ± 10°, corresponding to a transverse ve- 
locity of 171 km s""'^ ior D — 2 kpc. Their results are con- 
sistent with those reported by Wyckoff & Murray (1977) 
using historical images of the Crab Nebula. Though ~ 2 
kpc away, differential galactic rotation is negligible because 
the pulsar is within 5° of the Galactic anticenter direction. 
Caravco & Mignani pointed out that the proper motion 
and the symmetry axis of the inner X-ray nebula, based 
on ROSAT images, were essentially parallel. Arguing that 
the X-ray structure's symmetry is determined by the pul- 
sar's spin axis, they suggested that the (projected) spin 
axis and proper motion are parallel. The recent Chandra 
X-ray image (see Weisskopf 1999) confirms this picture 
which is also consistent with the optical work of Hester et 
al. (1995). While the morphology of the X-ray emission is 
less striictured and less symmetric than the Vela pulsar's 
nebula, there is a well-defined axis that indeed is parallel 
with the proper motion, to within the errors. 

Interpretation of the Crab pulsar's polarization is less 
clear than that of Vela, primarily because there is a multi- 
plicity of components whose relationship to the spin axis is 
ambiguous. However, radio polarization at 1.4 GHz (Mof- 
fett & Hankins 1999) has position angle (after correction 
for Faraday rotation) V'poi — —60° ± 10° for both the main 
pulse (MP) and interpulse (IP) components. At 5 GHz, 
the polarization angle for the MP is the same while the 
IP angle is shifted by 90° to 30° ± 10°. At 8 GHz, the 
MP is not visible while the IP has the same position angle 
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as at 5 GHz. Two new radio components (HFCl, HFC2) 
appear at 8 GHz that arrive later than the IP and have 
position angles ~ 85° ± 15°. At frequencies less than 1 
GHz, a 'precursor' component appears before the MP but 
its polarization angle is not available. At optical and UV 
wavelengths (Graham-Smith ct al 1998, 1996), the polar- 
ization angles at the ccntroids of the MP and IP are the 
same as in the radio (modulo 90° polarization-mode am- 
biguities). 

Empirically, the polarization vector at the ccntroid of 
the MP is parallel to the proper motion vector. The sim- 
plest interpretation is that the polarization mode in this 
pulse component is parallel to the projected spin axis. 
However, alternative models exist for the Crab pulsar's 
emission that place the emission region at sufficiently high 
altitude (in 'outer-gap' models) that rotational aberra- 
tion would break the coplanarity of the line of sight, spin 
and magnetic dipole vectors (e.g. Romani & Yadigaroglu 
1995). The geometry of the Crab's X-ray nebula may help 
clarify the emission geometry. 

2.3. The NS-NS Binary, B1913+16 

Analysis of the orbital elements, proper motion, and 
limits on the geodetic-precession angle leads to the con- 
clusion that a kick resulted from the most-recent super- 
nova that produced the observed pulsar's companion and 
constrains the angle between the kick and orbital plane 
(Wcx et al. 2000). The angle is < 5 — 10 degrees when 
the following hold: (1) all angular momenta (spin and or- 
bital) were aligned prior to the explosion (consistent with 
the canonical evolutionary history of double neutron star 
systems), (2) the kick timescale was much shorter than 
the orbital period (currently ^ 8 hr), and (3) the orbit 
did not evolve significantly since the explosion (see Fryer 
& Kalogera 1997). If the current determination of the 
hard-to-measure geodetic angle is ignored then the angle 
between kick and orbital plane is limited to < 30 degrees 
primarily by the condition that the binary remain bound. 

For all such analyses the following is important: the 
axis of the spin of the most-recently formed neutron star 
is unobserved, so the spin-kick orientation can be inferred 
only if one also assumes the neutron star spin was aligned 
with the progenitor's spin. Although plausible, there is no 
proof that the spin direction of the iron core of this neutron 
star progenitor matched the prcsupcrnova orbital angular 
momentum. We conclude that the results for B1913+16 
do not demand a different sort of kick mechanism than 
for the Crab and Vela pulsars. At the same time, if spin- 
kick alignment holds for B 19 13-1-16, if the angle between 
kick and orbital plane is < 30 degrees and if all angular 
momenta align, then the spin of the most-recently formed 
neutron star in that system is misaligned with the spin of 
the progenitor. 

2.4. Other Pulsars 

The alignments of proper motion vectors and projected 
spin axes for the Vela and Crab pulsars may be contrasted 
with what we know about other pulsars' geometries. For 
other pulsars, we have no information about the orienta- 
tion of the spin axis from nebular morphology. All that 
is known is based on radio polarization. For the 28 pul- 
sars (besides Vela) analyzed by Deshpande et al. (1999), 
there appears to be no systematic relationship between 



the polarization and proper motion vectors. We have used 
a likelihood analysis to compare different model relation- 
ships between these vectors. The null hypothesis, that 
they are independently and uniformly distributed on the 
sphere, was compared with two competing hypotheses: (1) 
that there is perfect correlation between the spin axis and 
either one polarization mode or its orthogonal counterpart; 
(2) a hybrid model where some objects have uncorrelated 
angles and others have perfectly correlated angles, modulo 
the 90 degree ambiguity. The likelihood analysis indicates, 
formally, that the best model is one where 80% of the ob- 
jects have uncorrelated angles while 20% have determin- 
istically related angles. However, the likelihood for this 
model is not significantly better than the null hypothesis. 
For the population of objects, we conclude, therefore, that 
there is no strong relationship between the polarization 
angle and the orientation of the proper motion. 

The lack of any empirically establishable correlation for 
a sample of pulsars should not be viewed as evidence that 
kicks and spins are usually unrelated. At least four effects 
can contaminate a potential correlation. (1) As discussed 
in §2.1, there are large uncertainties in using the polariza- 
tion angle to determine the spin axis; while the rotating 
vector model works well in fitting polarization angles ver- 
sus pulse phase in many cases, in others it does not. Also, 
aberration introduces an angular shift that depends on 
emission altitude, which may vary from object to object 
(Blaskiewicz, Cordcs & Wasserman 1991). "Orthogonal" 
polarization modes for some objects are not always pre- 
cisely orthogonal, raising the question that inference of 
the orientation of the spin axis from the polarization may 
be biased from the true value in some cases. About half 
the objects in Deshpande et al. 's sample are said to be 
affected by orthogonal mode flips. (2) Most neutron stars 
may be formed in binary systems and in such cases their 
present space velocities are a combination of disrupted or- 
bital motion and one or two kick velocities. Any spin- 
axis/kick velocity relationship c;an be diminished through 
this combination. While it is true that some pulsars have 
sufficiently large space velocities that any orbital contri- 
bution would be small, it is also true that other objects 
have sufficiently small space velocities that orbital motion 
may have contributed significantly. (3) Differential galac- 
tic rotation is important for distant objects and cannot be 
accounted for unless the distance is known accurately and 
the object has not moved far from its birth location. Seven 
objects in the Deshpande et al. sample have estimated dis- 
tances > 3 kpc and it is possible that others originated far 
enough away to be so affected. (4) The spin-velocity an- 
gle can be altered by acceleration in the Galactic potential 
for objects older than about 10 Myr. Five objects in the 
sample exceed this age. 

We stress that Deshpande et al. (1999) were careful to 
select objects whose polarization angle sweeps appeared 
to match the shape expected from the rotating-vector 
model. However, that model has sufficient parameters to 
fit any smooth position-angle curve even if that curve does 
not derive from the physical conditions that underly that 
model. A cogent example is, once again, the NS-NS bi- 
nary B1913+16. At 1.4 GHz, its pulse shape is dominated 
by two components that can be associated with a 'conal' 
beam; the position-angle curve indeed appears to follow 
a curve similar to that expected from the rotating-vector 
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model at the pulse phases of these components. How- 
ever, there is also contribution to the pulse from a 'core' 
component, as evidenced by its appearance at lower fre- 
quencies. In the 1.4 GHz polarization data, the position 
angle curve deviates markedly from the rotating-vector 
model at pulse phases corresponding to the core com- 
ponent (Cordes, Wasserman & Blaskiewicz 1990). This 
raises the suspicion that the polarization angle sweep of 
core components does not map the projection of the (as- 
sumed) dipolar field lines and thus does not allow deter- 
mination of the orientation of the spin axis. Given that 
core-component emission probably arises at or near the 
surface of the neutron star (e.g. Rankin 1990), it is plausi- 
ble that the magnetic field is distorted from a dipolar form 
there. Several objects in Deshpande et al. 's sample have 
profiles dominated by core components and several others 
show core components combined with conal emission. 

In summary, statistical studies do not allow one to con- 
clude that the angle between spin axes and kick directions 
is \miformly distributed on the sky. The Crab and Vela 
pulsars are not affected by these uncertainties because the 
Chandra X-ray data give an independent orientation vec- 
tor that is presumed to coincide with their spin axes. 

3. IMPLICATIONS FOR KICK MECHANISMS: NATAL KICKS 

Recent studies have focused on natal kicks imparted to 
the neutron star at birth. The strongest observational ev- 
idence that kicks must be natal comes from spin-orbit 
misalignment in the binary pulsar system PSR J0045- 
7319 (Kaspi et al. 1996) and PSR 1913-M6 (Kramer 1999) 
needed to produce the observed precessions. The alterna- 
tive, post-natal electromagnetic rocket mechanism (Harri- 
son & Tademaru 1975) will be discussed in §4. 

We discuss below the implications of aligned spin-kick 
in two classes of natal kick mechanisms (hydrodynamical 
driven and neutrino-magnetic field driven). 

3.1. Hydrodynamically Driven Kicks 

The first class relies on hydrodynamical perturbations 
in core collapse and supernova explosions. Numerical sim- 
ulations indicate that local hydrodynamical (convective) 
instabilities in the collapsed stellar core and its surround- 
ing shocked mantle (e.g., Herant et al. 1994; Burrows et 
al. 1995; Janka & Miiller 1994,1996; Keil et al. 1996; Mez- 
zacappa et al. 1998), which might in principle lead to 
asymmetric matter ejection and/or asymmetric neutrino 
emission, arc not adequate to account for kick velocities 
Z 100 km s"^ (Janka & Miiller 1994; Burrows & Hayes 
1996; Keil 1998). Global asymmetric perturbations of 
presupernova cores are required to produce the observed 
kicks hydrodynamically (Goldreich, Lai & Sahrling 1996, 
hereafter GLS: Burrows & Hayes 1996). GLS suggested 
that overstable g-mode oscillations in the presupernova 
core driven by shell nuclear burning may provide a nat- 
ural seed for the initial asymmetry. Calculations based on 
presupernova models of Weaver & Woosley (1993) indicate 

^Perturbations in the inner core (which collapses homologously) are damped during collapse (Lai 2000), thus the proto-neutron star is 
globally spherical. 

^The minimum rotation period of an isolated (i.e., decoupled from the envelope) core is about 1 s, and a real core must rotate signifi- 
cantly slower than this due to core-envelope coupling. Note that even if the newly formed neutron star rotates at the maximum breaJjup rate 
(Ps ^ 0.5 ms), the corresponding rotation period of the precollapse core is only ~ 11 s. 

^For example, the density, temperature structure of the rotating presupernova stellar model of Heger, Langer &; Woosley (2000), which does 
not include magnetic core-envelope coupling and thus overestimates the rotational effect, is not much different from the nonrotating model. 



that some g-modes are indeed overstable (Lai & Goldreich 
2000b: see Lai 1999), although uncertainties in the presu- 
pernova models preclude a definitive prediction of the per- 
turbation amplitudes. Alternatively, violent convection in 
the 0-Si burning shell may also produce asymmetric per- 
turbations (Bazan & Arnett 1998), although it is not clear 
whether such perturbations have sufhciently large scales. 
During collapse, the asymmetric perturbations seeded in 
the outer region of the iron core are amplified (by a fac- 
tor of 5-10) by gravity (Lai & Goldreich 2000a). ^ The 
enhanced asymmetric density perturbation leads to asym- 
metric shock propagation and breakout, which then gives 
rise to asymmetry in the explosion and a kick velocity to 
the neutron star (Burrows & Hayes 1996). 

Now let us consider how rotation might affect the kick. 
The low-order g-modes trapped in the presupernova core 
(M ~ 1.4M0, R ~ 1500 km) have periods of 1-2 seconds, 
much shorter than the rotation period of the core^, thus 
the g-modes are not affected by rotation. Also, since the 
rotational speed of the core is typically less than the speed 
of convective eddies (-1000-2000 km s'^, about 20% of 
the sound speed) in the burning shell surrounding the iron 
core, rotation should not significantly affect the shell con- 
vection either"^. (The convection zone provides an evanes- 
cent region for trapping the g-modes). Under these condi- 
tions, the development of large-scale presupernova (dipo- 
lar) asymmetry is not influenced by the core rotation. 

Even though the primary thrust to the neutron star 
(upon core collapse) does not depend on spin, the net kick 
will be affected by rotational averaging if the asymmetry 
pattern (near the shock breakout) rotates with the matter 
at period P shorter than the kick timescale Tkick- Let Vq 
be the kick velocity that the neutron star attains in the 
case of zero rotation, and let 9 be the angle between the 
primary asymmetry and the rotation axis. The expected 
components of kick along the rotation axis and perpendic- 
ular to it are (for Tkick ^ P) 



kick I 



= Vo cos 



ick_L 



V2P 
277 Tkick 



Vosinfil. (1) 



Thus the angle (3 between the kick vector Vkick and the 
spin vector ftg is given by tan/3 ~ 0.2(F/Tkick) tan6'. Typi- 
cally, the alignment between Vkick and fig will be achieved 
when Tkick » P- 

What is the kick timescale Tkick? In the standard 
paradigm of core-collapse supernovae, the shock wave gen- 
erated by core bounce starts at a radius of — 20 km (which 
encloses O.7M0), reaches 100 — 200 km in about 20 ms 
and stalls for hundreds of milliseconds until it is revived 
by neutrino heating (perhaps enhanced by vigorous con- 
vection in the mantle) and then moves outward again. As 
the shock breaks out asymmetrically, momentum is im- 
parted to the matter that will be incorporated into the 
neutron star on a timescale somewhat shorter than the 
delay. A reasonable estimate is Tkick ~ 100 ms which 



5 



since the cross section for z/g (z^e) absorption on neutrons 
(protons) depends on the local magnetic field strength, 
asymmetric neutrino emission can be produced if the field 
strengths at the two opposite poles of the star are differ- 
ent. To generate l^ick ^ 300 km using this mechanism 
would require that the difference in the field strengths at 
the two opposite stellar poles be at least 10^^ G (Lai & 
Qian 1998). 

A superstrong magnetic field may also play a dynam- 
ical role in the proto-neutron star. For example, it has 
been suggested that the magnetic field can induce "dark 
spots" (where the neutrino flux is lower than average) on 
the stellar surface by suppressing neutrino-driven convec- 
tion (Thompson & Duncan 1993). While it is difficult 
to quantify the kick velocity resulting from an asymmetric 
distribution of dark spots, a simple estimate indicates that 
a local magnetic field of at least 10^^ G is needed for this 
effect to be of importance. 

In this second class of mechanisms, the kick is imparted 
to the neutron star near the neutrinosphere^ on the neu- 
trino diffusion time, Tkick ~ 10 seconds. As long as the 
initial rotation period of the neutron star is much less 
than a few seconds, spin-kick alignment is naturally ex- 
pected [see eq. (1)]. Of course, while soft gamma re- 
peaters and anomalous X-ray pulsars may possess mag- 
netic fields in excess of 10^'' G ( "magnetars" ; e.g., Thomp- 
son & Duncan 1996; Vasisht & Gotthelf 1997; Kouve- 
liotou et al. 1998,1999), it is imclear (and perhaps \m- 
likely) that radio pulsars (and Vela and Crab in particu- 
lar) had B ^ 10^^ G at birth for these neutrino-magnetic 
field driven kicks to be relevant. 

3.3. Kick Induced Spin? 

For the two classes of kick mechanisms discussed in 
§3.1 and §3.2, spin-kick alignment requires that the proto- 
neutron star have a "primordial" rotation (i.e., with an- 
gular momentum coming from the presupernova core and 
possibly at a maximum rate if the kick is generated hy- 
drodynamically). How about the notion that the pulsar's 
spin is generated by off-centered kicks (Spruit & Phinney 
1998)? It is certainly true that even with zero precollapse 
angular momentum, some rotation can be produced in the 
proto-neutron star (Burrows et al. 1995 reported a rotation 
period of order a second generated by stochastic torques 
in their 2D simulations of supernova explosions): a kick 
Vicick = 3OOV300 km s~^, displaced by a distance s from 
the center, produces a spin Avg ^ 12V3oo(s/10 km) Hz, 
with the spin axis necessarily perpendicular to the kick 
direction (Burrows et al. 1995; Spruit & Phinney 1998). 
Aligned spin-kicks may be possible if the kick is the result 
of many small thrusts which are appropriately oriented 
(Spruit & Phinney 1998) — such a picture might apply if 
small-scale convection were responsible for the kick. How- 
ever, numerical simiilations indicate that such convection 
alone does not produce kicks of sufficient amplitude (Janka 
& Miiller 1994; Burrows & Hayes 1996; Keil 1998). As dis- 
cussed in §3.1 and §3.2, kicks of a few hundreds km s~^ can 

*In a prompt explosion, the shock wave exits the outer core (and the mass cut boundary) in about 20 ms; thus Tkick ~ 20 ms. Spin-kick 
alignment will result only when Ps <S 0.2r^^ ms, where r2 can be as small as 0.5, corresponding to the radius of the hot proto-neutron star. 

^Here we do not discuss kick mechanisms based on nonstandard neutrino physics [e.g., Kusenko & Segre (1996); Akhmedov et al. (1997); 
but see Janka & Raifelt (1998) for critiques on these papers]. 

®The neutrinosphere of the proto-neutron star decreases from ^ 50 km to near its final value of 10-20 km in the first 100 ms after bounce; 
see Burrows & Lattimer 1986. 



is the shock travel time at speed of 10 km s~^ across 
^ 1000 km, the radius of the mass cut enclosing 1.4A'/q 
(e.g., Bethe 1997; Burrows & Hayes 1996). For spin-kick 
alignment, the rotational period of the perturbation pat- 
tern at r = 100r2 km must be much less than Tkick- This 
corresponds to a final neutron star (i? = 10 km) spin pe- 
riod Ps <C Tkick/(100 r2) ^ 1 ''2^^ ms. Note that since the 
shock wave travels from ~ 100 km to ^ 1000 km while 
the kick is being imparted to the star (e.g.. Burrows & 
Hayes 1996), ^2 = 1 corresponds to the maximum spin 
period of the final neutron star. We thus conclude that 
if rotation is dynamically unimportant for the core col- 
lapse and explosion (corresponding to Pg ^ 1 ms), then 
rotational averaging is inefficient and the hydrodynamical 
mechanism does not produce spin-kick alignment.* 

The discussion above is based on the standard picture 
of core-collapse supernovae, which is valid as long as ro- 
tation does not play a dynamically important role (other 
than rotational averaging) in the supernova. If, on the 
other hand, rotation is dynamically important, the ba- 
sic collapse and explosion may be qualitatively different 
(e.g., core bounce may occur at subnuclear density, the 
explosion is weaker and takes the form of two-sided jets; 
Monchmeyer et al. 1991; Rampp, Miiller & Ruffert 1998; 
Khokhlov et al. 1999; Fryer & Heger 1999). The possibility 
of a kick in such systems has not been studied, but it is con- 
ceivable that an asymmetric dipolar perturbation may be 
coupled to rotation, thus producing spin-kick alignment. 

3.2. Neutrino - Magnetic Field Driven Kicks 

The second class of kick mechanisms relies on asymmet- 
ric neutrino emission induced by strong magnetic fields. 
^ The fractional asymmetry e in the radiated neu- 
trino energy required to generate a kick velocity Vkick is 
e = MVi,ickc/Etot (= 0.028 for Hick = 1000 km s-\ 
M = 1.4 A/q and total neutrino energy radiated Eioi = 
3 X lO^'' erg). In the magnetized nuclear medium, the neu- 
trino scattering/absorption opacities depend asymmetri- 
cally on the directions of neutrino momenta with respect 
to the magnetic field axis. (This is a manifestation of par- 
ity violation in weak interactions.) Asymmetric neutrino 
flux can be generated in the outer region of the proto- 
neutron star (i.e., above the neutrino-matter decoupling 
layer, but below the neutrinosphere) where the neutrino 
distribution deviates from thermal equilibrium (Vilenkin 
1995; Arras & Lai 1999a,b). Arras & Lai (1999a,b) found 
that averaging over all neutrino species, the total asymme- 
try in neutrino flux is e ~ 0.1 B-^^E'^ + O.OO2B15/T (the 
first term comes from the effect of quantized states of elec- 
trons in the absorption opacity, and the second term comes 
from nucleon polarization), where B = 10"'^^i?i5 G is the 
magnetic field strength, E^, is the mean energy (in MeV) of 
VgS and T is the temperature (in MeV) in the decoupling 
layer. The resulting kick velocity is \4ick ^ 50 Sis km s^^. 
Kicks of a few hundreds km s^^ would require the proto- 
neutron star to possess an ordered component of magnetic 
field with magnitude greater than 10^'"^ G. Alternatively, 
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be generated if global dipolar symmetry is broken (either 
due to hydrodynamical perturbations or due to magnetic 
fields). 

4. POST-NATAL KICK: ELECTROMAGNETIC ROCKET 
EFFECT 

Harrison & Tademaru (1975) show that electromagnetic 
(EM) radiation from an off-centered rotating magnetic 
dipole imparts a kick to the pulsar along its spin axis. 
The kick is attained on the initial spindown timescale of 
the pulsar (i.e., this really is a gradual acceleration), and 
comes at the expense of the spin kinetic energy. We have 
reexamined this effect and found that the force on the 
pulsar due to asymmetric EM radiation is larger than the 
original Harrison & Tademaru expression by a factor of 
four. If the dipole is displaced by a distance s from the 
rotation axis, and has components ^ip,^^,fiz (in cylindri- 
cal coordinates), the force is given by (to leading order in 
ns/c) 



^-15'- 



(2) 



(The sign is such that negative F implies Vkick parallel 
to the spin O.) The dominant terms for the spindown 
luminosity give 



3^ 



2a 



5c2 



(3) 



For a "typical" situation, fip ^ fiip ^ the asymme- 
try parameter e = F/{L/c) is of order OA{ils/c). For a 
given n, the maximum Cmax = \/0.4 — 0.63 is achieved 
for fJ.p/fiz — and fi^/fJ-z — VOA{Qs/c). From AIV — 
e{L/c) = — e(/ilil)/c, we obtain the kick velocity 



ykick^ 260i?^o( — 



oTi 



IkHz 



1- - 



km s 



km s (4) 



where R — 10i?io km is the neutron star radius, Vi is 
the initial spin frequency, v is the current spin frequency 
of the pulsar, and e = (O? — Jedi^^; in the sec- 

ond equality, we have adopted the "typical" situation 
{fip = = ^z) so that e = OA{Qs/c). For the "opti- 
mal" condition, with /ip = 0, /i^Z/iz — ^0.4 (Ois/c), and 
e = Vol [2n,n/{n^ + n^)] , we find 



^ki^r^ - 1400 



10 



VI kHz/ 



4.66 ( — 



-1 ^ 
tan — 



km s 



(5) 



{v ^ 100 Hz) neutron star can potentially lose signif- 
icant angular momentum to gravitational waves gener- 
ated by unstable r-mode oscillations (e.g., Andersson 1998; 
Lindblom, Owen & Morsink 1998; Owen et al. 1998; An- 
dersson, Kokkotas & Schutz 1999; Ho & Lai 2000). If 
gravitational radiation carries away the rotational energy 
of the neutron star faster than the EM radiation does, 
then the electromagnetic rocket effect will be much di- 
minished^. In the linear regime, the r-mode amplitude 
a ^/R (where ^ is the surface Lagrangian displacement; 
see the references cited above for more precise definition 
of a) grows due to gravitational radiation reaction on a 
timescale tgrow — 19 (j^/l kHz)~^ s. Starting from an ini- 
tial amplitude ^ 1, the mode grows to a saturation 
level ttsat in time tgrow ln(asat/Q!i) during which very little 
rotational energy is lost. After saturation, the neutron star 
spins down due to gravitational radiation on a timescale 



TGR 



100 a: 



GR 



(ikHz) 



(6) 



(Owen et al. 1998). By contrast, the spindown time due 
to EM radiation alone is 



tem 



10^ B7^^ 



EM 



(ikHz) 



(7) 



where B13 is the surface dipole magnetic field in units of 
10^'^ G. Including gravitational radiation, the kick velocity 
becomes 



kick 



Mc 



TGR 



'''GR + '''EM 



--'o(n)(TSi) 



In 

[3 



1 



km s 



(8) 



where in the second equality we have replaced e by con- 
stant mean value e, and is defined by 







Vtgr 



( "sat \V l^r Y 

V 10-2-5 J V IkHz J ^1 



(9) 



Thus if the neutron star was born rotating at ^ 1 kHz, 
it is possible, in principle, to generate spin-aligned kick of 
a few hundreds km or even 1000 km s~^. 

Equations ^ and (H) assume that the rotational en- 
ergy of the pulsar entirely goes to electromagnetic radi- 
ation. Recent work has shown that a rapidly rotating 

^Gravitational radiation can also carry away linear momentum, but 



For /3 <C 1, equation (||) becomes for P ^ I, the kick 
is reduced by a factor 1//3. 

Clearly, for the EM rocket to be viable as a kick mech- 
anism at all requires P 1. The value of asat is unknown. 
Analogy with secularly unstable bar-mode in a Maclau- 
rin spheroid implies that asat '^-^ 1 is possible. It has been 
suggested that turbulent dissipation in the boundary layer 
near the crust (if it exists early in the neutron star's his- 
tory) may limit asat to a small value of order lO-^-iQ--^ 
(Wu, Matzner & Arras 2000). The theoretical situation is 
not clear at this point. 

The EM rocket effect has not been popular in recent 
years because empirical tests using radio polarization and 
proper motion have not yielded a positive correlation and 
because it is widely thought that radio pulsars are born 
with spin periods of 0.02-0.5 s (e.g., Lorimer et al. 1993). 
However, as discussed in §2.3, the polarization-proper mo- 
tion correlation might not be detectable even if kicks are 

the effect for a neutron star is negligible. 
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always aligned with spin axes. The recent discovery of the 
16 ms X-ray pulsar (PSR J0537-6910) associated with the 
Crab-like supernova remnant N157B (Marshall et al. 1998) 
implies that at least some neutron stars are born with spin 
periods in the millisecond range. In particular, a recent 
analysis of the energetics of the Crab Nebula suggests an 
initial spin period ~ 3 — 5 ms followed by fast spindown on 
a time scale of 30 yr (Atoyan 1999). As for the Vela pulsar, 
the energetics of the remnant do not yield an unambiguous 
constraint on the initial spin (S. Reynolds, private com- 
munication). Thus it seems prudent to consider the EM 
rocket effect as a possible kick mechanism. We do note, 
however, that such slow post-natal kick may have difficulty 
explaining some properties of binary pulsar systems, such 
as the spin-orbit misalignment in PSR J0045-7319 (Kaspi 
et al. 1996) and PSR 1913-hl6 (Kramer 1999) needed to 
produce the observed precessions. ^ Similarly, a slow kick 
(with Tkick ^ ^orb) may be inconsistent with the neutron 
star binary populations (e.g., Dewey & Cordes 1987; Fryer 
& Kalogera 1997; Fryer et al. 1998). Some of these issues 
will be addressed in a future paper. 

5. CONCLUSION 

Motivated by the apparent alignment between the spin 
axis and the proper motion for the Vela and Crab pul- 
sars as revealed by recent Chandra observations and by 
analysis of radio polarization (§2), we have examined in 
this paper the questions of how rotation affects kick and 
whether spin-kick alignment is generally expected in dif- 
ferent classes of theoretical models of neutron star kicks. 
We find that 

(1) For hydrodynamically driven kicks, spin-kick align- 
ment is possible only if rotation plays a dynamically impor- 
tant role in the core collapse and explosion, corresponding 



to neutron star initial spin period P ^ I ms (§3.1); 

(2) For neutrino-magnetic field driven kicks, spin-kick 
alignment is almost always produced (§3.2); but for these 
class of mechanisms to be relevant, the proto-neutron star 
must have magnetic field stronger than 10"'^^ G; 

(3) Post-natal kicks due to the electromagnetic rocket 
effect always lead to spin-kick alignment; even with our 
upward revised (by a factor of 4) force expression com- 
pared to the original Harrison- Tademaru formula, for this 
mechanism to be of interest, the proto-neutron star must 
rotate at frequency ^ 1 kHz, and gravitational radiation 
(via r-mode instability) driven spindown must be less ef- 
ficient than magnetic braking — this puts a constraint 
on the magnetic field and r-mode amplitude in the young 
neutron star (§4). 

Regardless of which kick mechanism is at work, if spin- 
kick alignment turns out to be a generic feature, then neu- 
tron stars must have non-negligible (and perhaps maxi- 
mum) initial rotation. Theoretically, it has been argued 
that efficient core-envelope coupling via magnetic stress 
in the presupernova star results in negligible "primordial" 
angular momentum for the neutron star (Spruit & Phin- 
ney 1998), although there are different views (e.g., Livio & 
Pringle 1998). If, on the other hand, future observations 
show that spin-kick alignment is not a generic feature, then 
we must conclude that the alignment in Vela and Crab pul- 
sars is a coincidence, and the origin of neutron star kicks 
lies in hydrodynamics. 

This work is supported in part by NASA grants NAG 
5-8356 and NAG 5-8484, by NSF grants AST 9819931 and 
AST 9986740, and by a research fellowship (to D.L.) from 
the Alfred P. Sloan foundation. 
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